Engineering spin-waves in a high-spin ultracold Fermi gas 
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We report on the detailed study of multi-component spin- waves in an s = 3/2 Fermi gas where 
the high spin leads to novel tensorial degrees of freedom compared to s= 1/2 systems. The exci- 
tations of a spin-nematic state are investigated from the linear to the nonlinear regime, where the 
tensorial character is particularly pronounced. By tuning the initial state we engineer the tensorial 
spin- wave character, such that the magnitude and sign of the counterflow spin- currents are effec- 
tively controlled. A comparison of our data with numerical and analytical results shows excellent 
agreement. 



PACS numbers: 05.30Fk, 03.75Ss, 75.30Ds 

Spin-interaction driven phenomena are crucial for the 
behavior of many quantum systems, e.g., ferromagnets 
[1] and high-temperature superconductors and are 
also relevant in spintronics applications [3]. Apart from 
condensed matter systems with an electronic spin of 
5 = 1/2, dilute atomic gases show a wealth of novel spin- 
excitations, where the spin is provided by the internal 
hyperfine structure of the atoms. Pioneering experi- 
ments with hydrogen ^ and helium [5 showed the ex- 
istence of transversal spin- waves, which arise from in- 
trinsic spin-exchange interactions [BHS]- Longitudinal 
spin- waves in two-component mixtures have been demon- 
strated in non-condensed bosonic ^^Rb [9 . For weakly in- 
teracting fermions, slow spin-currents were reported near 
the zero-crossing of a Feshbach-resonance [10] and the 
interaction-induced damping of dipole oscillations was 
studied |TT]. Prominent examples for strongly interact- 
ing fermionic systems are the BEC-BCS crossover [I2j , 
miscibility of spin mixtures [13], and the quest for itin- 
erant ferromagnetism [MHTT] . In contrast to conven- 
tional two-component systems, the hyperfine structure 
of many atoms also allows for spinor gases with 5 > 1/2, 
which offer a whole new set of possibilities to study spin- 
dependent phenomena ^1 . This includes the exis- 
tence of spin-nematic states [20l [21^ , spin-changing colli- 
sions [221 [23], spontaneous domain formation [24], SU(N) 
degenerate ground states [2F-29^, and novel superfluid 
phases |30l|3l]. Especially for fermionic atoms, s = 3/2 
constitutes the simplest realization of a high-spin system 
and is thoroughly studied theoretically, being a model 
system for all higher spins [32 , 33 . 

In this Letter, we investigate tunable spin- waves in 
a Fermi gas with pseudo-spin s = 3/2 in the quantum 
degenerate regime, which constitute fundamental collec- 
tive excitations in these systems. Spin-waves in such 



high-spin systems are predicted to exhibit very complex 
and novel properties (see Refs. [34l [35] for bosonic sys- 
tems): Higher spins require additional tensor compo- 
nents in addition to the well-known spin-vector used to 
describe s = l/2 system, such as nematic and octupole 
operators in the case of s = 3/2. The coupling between 
these tensor components strongly influences the spin- 
wave properties and leads to completely new phenom- 
ena. We have investigated spin-wave excitations ranging 
from the linear regime, where only spatial dipole exci- 
tations are present to the non-linear regime, where the 
nematic spin-component predominantly couples to the 
quadrupolar spatial mode. The change of the spin- wave 
composition is accompanied by a change of the oscillation 
frequency, which has a minimum in the linear regime. In 
a second set of experiments, we engineer the initial state 
coherences. By this, the spin-current for two of the four 
components can be reversed changing the spin- wave char- 
acter from spin-octupole to spin- vector. All results are in 
excellent agreement with numerical calculations based on 
a semiclassical mean-field theory (SMFT; a collisionless 
Boltzmann equation with full spin coherence) [35H38] . 
Our results demonstrate the high degree of control that 
can be exerted on spin-waves in high-spin Fermi gases 
and thus pave the way towards novel schemes for atom 
spintronics using the intrinsic high spin. 

Our measurements are performed in a quantum de- 
generate gas of ^^K in the / = 9/2 hyperfine manifold. 
We initially evaporate a balanced mixture of |m = l/2) 
and |m = — 1/2) to quantum degeneracy in an optical 
dipole trap [36]. The final trapping frequencies are 
^x,y,z = 27r X (70, 70, 12) Hz. At this point, we apply a 
radio- frequency (rf) pulse to create a coherent superpo- 
sition with the states |±3/2) [Fig.[T](a)]. We initialize the 
spin- waves by applying a small magnetic field gradient up 
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FIG. 1. (a) Single-particle density matrix for an incoherent 
superposition of |l/2) and | — 1/2) (left) and the resulting co- 
herent superposition of all four components after a resonant 
rf-pulse (right). Diagonal elements Wa are real and repre- 
sent the populations. Off-diagonal elements are complex num- 
bers Wij = |Wij|e^^^j and include the phase ^ij between differ- 
ent components. Plotted is only the absolute value |Wij| (b) 
Sketch of the local phase across the Fermi gas after pulses with 
different magnetic field gradients, (c) m = component of the 
Z = 0, 1, 2, 3 tensor operators T™ for s = 3/2 in comparison to 
the corresponding spherical harmonics y^™. 



to a few G/m for 10ms, which leads to a phase spiral for 
coherent superpositions of different spin-components as 
sketched in Fig. [ijb). While these coherent superposi- 
tions are initially still spin-polarized locally, the phase- 
twist allows for interactions in a trapped gas where the 
external potential induces spatial dynamics [37]. In gen- 
eral, the resulting mean-field interaction couples the spin 
degrees of freedom to different modes of the external trap 
leading to the emergence of spin- waves. We detect the 
spin-current using absorption imaging either in situ or 
after 18.5 ms time-of-flight (TOF) with a Stern-Gerlach 
separation of the spin-components [36]. In Fig. |2|a) we 
show a typical example for a 5 = 3/2 spin- wave initial- 
ized by a 10ms gradient of AB = 3AG/m. The mea- 
surements reveal oscillatory spin-currents in all four spin- 
components. We observe a time-independent total den- 
sity, meaning that these spin-waves constitute counter- 
flow spin-currents without an accompanying net mass 
transport. In particular, note the inverted flow direc- 
tion of |l/2) (1-1/2)) with respect to the |3/2) (|-3/2)), 
which is a clear indication of the new tensorial degrees of 
freedom as discussed later. 

Let us at this point introduce a description of spin in 
the language of tensors which transform invariant under 
rotations, to obtain a deeper insight into the underlying 
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FIG. 2. (a) In situ time evolution of all four spin- 
components after a 10 ms pulse with a magnetic field gradient 
of Ai? = 3.4G/m. Shown are the colum densities at different 
times after the excitation, (b) Deviation from the initial pop- 
ulation of the m = component of the vector (/ = 1), nematic 
(/ = 2) and octupole (/ = 3) component. The vector and oc- 
tupole component show spatial dipole oscillations, while the 
nematic component clearly exhibits quadrupolar spatial os- 
cillations. (c,d) Numerical calculation for the parameters of 
(a,b). 



physical processes. Most common are the spherical har- 
monics for orbital angular momentum [Fig.jljc)] and the 
spin- vector S oc {ax^ az) where are the Pauli ma- 
trices. The latter are conveniently used as a basis set 
to describe any (pseudo) 5 = 1/2 system. To describe 
the physics of larger spins, however, it is necessary to 
include higher order tensors. In a spin 3/2 system, as 
discussed here, in addition to the spin- vector (/ = !), the 
spin-nematic tensor (/ = 2) and the spin-octupole tensor 
(/ = 3) must be included (see Fig. [TJc) and [36]). In par- 
ticular, this basis set automatically fulfills all possible 
rotation symmetries which simplify the equations of mo- 
tion. The initial state in Fig. [2] is explicitly chosen to be 
purely nematic to exploit these symmetries in the exper- 
iments. The time evolution of the / = 1,2,3 tensors for 
the measurements of Fig.[2|a) is shown in Fig.[2|b). Note 
the predominantly quadrupolar spatial dynamics of the 
spin-nematic component, qualitatively different from the 
dipolar oscillations in the spin-vector and spin-octupole 
component. This results from a self-decoupling of the 
nematic component due to the rotational symmetry and 
our properly chosen initial state, as explained later. 

High-spin spin-waves can be theoretically treated 
within SMFT j34H36] . In the tensor notation, this is 



a generalization of the usual spin-vector approach for 
s = 1/2 systems [lOl [37H4Q] , which allows to understand 
the intriguing observations from Fig. [2] The system is 
described by a single-particle Wigner function W{z^p) 
in one spatial dimension and the equation of motion 
takes the form of a Boltzmann-equation in the collision- 
less regime. To zeroth order in the SMFT one gets for a 
s = 3/2 gas in a spin-independent external harmonic trap 

dtW{z,p) = doW{z,p) + ^[W{z,p), V{z)] . (1) 

Here, 14in(^)= /Z)ki(^kinm - U]^mni)W]^i{z,p)dp is the 
effective mean-field interaction with the spin-dependent 
coupling constants [/yki [36^, do = {—p/mdz -\- mujlzdp) 
captures the time evolution due to the harmonic trap 
and kinetic energy, m is the mass of ^^K, and [., .] in- 
dicates the commutator in spin space. In the simula- 
tions, also the first order gradient expansion in the in- 
teractions are taken into account, leading to very small 
deviations only [36 . For a s = 3/2 system, the Wigner 
function in spin space is a 4 x 4 single-particle density 
matrix (SPDM), where the diagonal elements Wu rep- 
resent the absolute population of the spin-components 
and the off-diagonal elements Wij = | Wij |e^^^j represent 
the single-particle-coherences between different compo- 
nents. To induce a time evolution of the populations 
VFii, it is sufficient to spatially vary the phases ^ij of the 
off-diagonal elements , since they are coupled via the 
commutator in ([T]). 

The equation of motion ([T]) can be cast into the tensor 
notation introduced above by decomposing the Wigner 
function (mean-field potential) into tensor components 
Wi (Vi) [41,: 

dtWo = doWo , 

in 

dtW2 = doW2 + 4 {[W2, V, + Vs] + [Wi + W3, V2]) , 

dtWs = doWs + 4 ([W3, Fi] + [Wi + W3, V3] + [W2, V2]) . 

(2) 

The structure of Q leads to several important conse- 
quences. First, the total density Wq is not affected by the 
phase spiral, since it contains no off-diagonal elements. 
Therefore, the total density remains constant as we ob- 
serve in the experiment. Second, the dynamics of the 
nematic tensor W2 is proportional to the vector and oc- 
tupole components {Wi and IV3), while the latter two 
have terms proportional solely to the nematic / = 2 com- 
ponent. A purely nematic initial state as in Fig. |2] does 
therefore not support nematic excitations to first order. 
In the highly nonlinear regime, where vector and octupole 
excitations possess a large amplitude, nematic excitations 
are created via higher-order processes. This leads to the 
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FIG. 3. (a) Frequency and (b) oscillation amplitude of 
spin-waves excited with different magnetic field gradients at 
cjz = 12 Hz. Negative amplitudes denote an inverted initial 
direction. Solid lines are numerical data for each component. 
All error bars solely correspond to fit errors, representing one 
standard deviation. The experimental amplitudes which are 
taken after TOF and the numerical amplitudes calculated in 
situ are rescaled onto each other by a global factor, (c) Sketch 
of the phase windings across the atom cloud for different gra- 
dient pulses. 



spatially quadrupolar dynamics of the nematic compo- 
nent, which is qualitatively different from the spatially 
dipolar oscillations of the vector and octupole compo- 
nent [Fig. [2|b)]. Figure [2|c,d) show numerical results 
using ([1]) for the exact experimental parameters. Note, 
that these are in excellent agreement with the experimen- 
tal results demonstrating the capability to understand 
complex spin- waves via the demonstrated tensorial spin- 
decomposition. 

To analyze the behavior of the system for different exci- 
tation amplitudes, we applied different gradient strengths 
during the initialization of the spin- wave [42] . This corre- 
sponds to a change of the initial phases in the SPDM, 
while the initial coherence amplitudes \ Wij \ are kept con- 
stant. In Fig. [3] experimental results are compared to nu- 
merical calculations and show very good agreement: For 
small gradients, the frequency is amplitude-independent 
and the amplitude rises approximately linear with the 
gradient strength. For large gradients, the frequency ap- 
proaches the trapping frequency and is again only weakly 
depending on the excitation amplitude. For intermediate 
gradients, the system behaves strongly nonlinear which 
results in an amplitude-dependent oscillation frequency. 
In the regime of small gradients, corresponding to small 
excitation strengths one can apply the method of mo- 
ments [43 up to linear order in z and p, which describes 
pure spatial dipole oscillations ^ . Their oscillation fre- 
quency for the present initial state can be derived to be 
- cjmf 7 where cjmf is the mean-field inter- 



mf 



action energy as defined in Ref. [36] . This frequency is de- 
termined by a competition between the trap induced spa- 
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FIG. 4. (a) Spatial oscillation amplitude of the spin- wave ex- 
citations for different initial coherences but equal populations 
of the four components at AB = 3.6 G/m. Solid lines show the 
initial spin-wave amplitude extracted from numerical calcula- 
tions, (b) Exemplary SPDMs for different initial coherences, 
(c) Amplitude of the dipole and octupole tensor components 
[36]. (d) Vector and octupole tensors evaluated in (c). All 
error bars solely correspond to fit errors, representing one 
standard deviation. The experimental amplitudes which are 
taken after TOF and the numerical amplitudes calculated in 
situ are rescaled onto each other by a global factor. 



tial oscillations and the mean-field rotation of the spin. 
For our parameters, we calculate = 27r x 2.3 Hz in good 
agreement with the experimental results. The linearized 
equations of motion also confirm the pure vector and oc- 
tupole character of the dipole excitations for a perfectly 
nematic initial state as discussed above. 

For large gradients, the frequency slowly approaches 
the trapping frequency: In this regime, the phase spi- 
ral only leads to an initial strong kick but is averaged 
out dynamically on timescales 2i\ ju^ |37j. Consequently, 
the mean-field potential no longer affects the subsequent 
oscillation and the corresponding frequency approaches 
the trap frequency with increasing winding number of 
the phase spiral, equivalent to an increasing gradient. 

All measurements discussed so far were performed with 
the same purely nematic initial state. By modifying the 
rf-pulse sequences used for the preparation of the ini- 
tial state we can control the amplitude of the coherences 
\Wij I and populations VKn in the SPDM. By this the mul- 
tipole decomposition of the initial state can be widely 
tuned and allows for the initialization of e.g. pure vector 
or nematic initial states, which in turn results in differ- 
ent spin and spatial characteristics of the emerging spin- 



wave. Following this direction, we performed a second 
set of experiments, where we engineered the spin- wave 
excitations by keeping the population of all four spin- 
components constant but changing the initial coherence 
amplitudes ll^ijl [Fig. Qb)]. Note, that this is comple- 
mentary to the results shown in Fig.jsj where we changed 
the phase 6>ij of the coherence by using different gradi- 
ent strengths. Figure [4]^ a) shows the resulting oscilla- 
tion amplitude for all four spin-components depending on 
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At c~0.5 the system changes 



its qualitative behavior where the |±l/2) components in- 
terchange their oscillation direction. Using the tensor 
notation, the spin-wave at small c is dominated by the 
spin-octupole, where neighboring spin-components have 
an inversed sign and therefore oscillate in opposed di- 
rections. At large c, the oscillation becomes spin-vector 
dominated, where spin-components with the same sign of 
magnetization oscillate in the same direction [Fig. |4]^c)]. 
The anew increase of the octupole amplitude at large 
c is due to higher order spatial excitations, possible in 
the nonlinear regime, where the measurements were per- 
formed. At c~0.5, the vector and octupole component 
contributions mutually cancel, leading to a vanishing 
spin-current in the |±l/2) components. Again the nu- 
merical calculations describe the engineered spin-waves 
very well. 



In conclusion, we have thoroughly investigated the 
physics of tunable spin- waves in a high-spin Fermi gas. 
We have generalized the spin- vector description of s = 1 /2 
systems by introducing higher order spin-tensors and 
used the corresponding decomposition of the equations of 
motion to explain the novel emerging spin-wave charac- 
teristics in a high-spin system. In particular, we observe a 
self-decoupling of the spin- nematic component, suppress- 
ing nematic spin- waves. Further, we have analyzed the 
spin-wave excitation spectrum for different phase spirals, 
corresponding to different excitation strengths ranging 
from the linear deep into the nonlinear regime. Finally, 
we demonstrated how to control the multipole charac- 
ter of spin-waves which leads to a reversal of the re- 
sulting counterflow spin-current of two spin-components. 
Our results constitute the first study of collective spin- 
spatial excitations in high-spin bulk Fermi gases and 
demonstrate the controlled manipulation of atomic spin- 
currents. This paves the way towards novel schemes for 
spintronics in ultracold atomic gases, using the intrinsic 
high spin as a valuable resource. 
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SUPPLEMENTAL INFORMATION 

In this supplemental material we discuss the 
preparation of our fermionic quantum gas (SI), 
the detection after TOF and in situ (S2) and the 
data analysis (S3). The theoretical model is in- 
troduced in (S4) and the linearized calculation is 
presented in (S5). For the tensor expansion in 
the 5 = 3/2 system see (S6). 

SI. PREPARATION OF FERMIONIC QUANTUM 
GASES 

We sympathetically cool about N = 2 x 10^ spin- 
polarized ^^K atoms in the state / = 9/2 and |m = 9/2) 
to 0.1 Tp using ^''Rb in a magnetic trap. Afterwards 
we transfer the atoms to a crossed circular-elliptical op- 
tical dipole trap with A = 812nm. The radii are 
'^x,y = 120 jam for the circular beam and = 70 /im and 
Wz = 280 jam of the elliptical beam, where the tight focus 
is in the vertical direction. Using rf-pulses and rf-sweeps, 
an equal mixture of two the hyperfine states |l/2) and 
| — 1/2) is prepared and e vapor at ively cooled in the trap 
by a 2s exponential intensity ramp. The final trapping 
frequencies are typically 0Jx,y,z = 27r x (70, 70, 12) Hz with 
AT = 3.5 X 10^ particles at 0.25 Tp. 

S2. DETECTION AND ANALYSIS OF SPIN 
COMPONENTS 

To obtain the experimental data, we used two differ- 
ent detection protocols, which have complementary ad- 
vantages. In most of the experiments, we applied time of 
fiight imaging, where we switch off all optical potentials, 
and perform a Stern-Gerlach separation of the different 
spin states within the free expansion time of 18.5 ms. The 
atoms are detected via resonant absorption imaging and 
the center-of-mass is calculated for each of the separated 
clouds individually. This method has the advantage that 
all spin-components can be detected simultaneously, but 
the disadvantage that all spatial modes but the dipole 
mode are effectively washed out, since the TOF mixes 
momentum and spatial components of these modes. 

To obtain more information on the spatial modes of 
the excited spin-waves we employed an in situ detection 
protocol, allowing us to observe, e.g., the quadrupole os- 
cillations in Fig. 2(b). Instead of separating the differ- 
ent spins by a Stern-Gerlach procedure in TOF, we use 
microwave pulses at 1.3 GHz, to transfer all but one sin- 
gle component to the f = 7/2 manifold of ^^K, which is 
off-resonant to the detection light, and consequently the 
transferred atoms do not appear in the absorption im- 
ages. Afterwards we switch off all magnetic fields and 
optical potentials for 1 ms to detect the atomic sample. 



To record the time evolution of all four components, it is 
therefore necessary to repeat the full measurement four 
times. 

S3. ANALYSIS OF EXPERIMENTAL DATA 

For every time step of the spin- wave dynamics, we de- 
termine the center-of-mass of each of the spin compo- 
nents. We extract the oscillation frequencies uj by fitting 
an exponentially damped cosine of the form 

COM(t) = A exp(-n) cos(cjt + c^) + C , (3) 

with oscillation amplitude A, damping rate F, a phase 
shift $ and a constant offset C. For the dipole and oc- 
tupole fraction in Fig. 4(c), we use a slightly different 
method. We take the spatial average of the absolute value 
of the projected multipoles M for each time step as 

^W = ^El^(^'^)l (4) 

where p is the pixel index of the image and A^p is the num- 
ber of pixels. By this procedure, the sign-informations of 
the projection are lost and therefore, we fit the resulting 
time series by 

M{t) = A exp(-rt) cos{ujt + <l>)^ + C (5) 

and report the amplitude A. 

S4. KINETIC EQUATION FOR A TRAPPED 
HIGH SPIN FERMI GAS 

We describe the multi-component fermionic system us- 
ing the single-particle Wigner-function, with its time evo- 
lution given by the semiclassical kinetic equation 

dtW{z,p) = doW{z,p) + ^ [W{z,p),V{z)] 

+ ^{dpW{z,p),d,V{z)} (6) 

with the short notation dn = — —dz -\- muj'^zdr} for the 
single-particle part of motion. We do not explicitly in- 
clude the magnetic field gradient, since it is used for the 
creation of the spin- wave but not necessary for its prop- 
agation once created. The derivation of this equation is 
outlined in previous works [371 131] • One basically starts 
from the full mean-field dynamics of the Wigner function 
and applies a semiclassical approximation by neglecting 
terms involving higher order derivatives with respect to 
position and momentum. To leading (zeroth) order, the 
mean-field interaction gives rise to the commutator and 
the next (first) order is given by the anticommutator. 
We have neglected here even higher orders terms. The 
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important term here is the commutator which is absent 
in the spinless case and describes interaction driven co- 
herent spin dynamics. It is also dominant with respect 
to the anticommutator and we neglect the latter in our 
analytical studies (while we keep it in the numerical sim- 
ulations). The kinetic equation (|6| is nonlinear, since 
the mean-field potential depends on the Wigner-function 
itself and its matrix elements are defined as 



Vmn{z) = / dp^{Uklr, 



Ukmnl)Wkl{z,p) (7) 



with the coupling constants defined as 
Uijki = El=o\... 9s Ei=-siSM\ik){SM\jl) for 
arbitrary spin s with the corresponding Clebsch-Gordan 
coefficients {SM\ik) = {SM\s^i^s,k). Here gs denotes 
the interaction strength for the scattering of two par- 
ticles with total spin S. For a real spin-3/2 system 
there are only scattering channels for 6' = 0,2, but in 
the numerical simulations of (|6| we take into account 
the values of Uijki for spin 9/2 in the sub-manifold of 
i,j,k,l= ±3/2, ±1/2, which also depend on ^4,6,8 [23] . 
Since in the experiment spin-waves along the z-axis are 
created, an effective one-dimensional kinetic equation is 
sufficient to describe these spin- waves. This is obtained 
by integrating out the transversal degrees of freedom 
from the three-dimensional case. The coupling constants 
for the 3D case, g'^ = ^^^as with s-wave scattering 
lengths a^, are modified by the local transversal density 

profile as gs = (/dJfXn[x^^^^^^^^ >< ^s- This modification 
depends on z but only weakly so we approximate it with 
the central density n(x,?/,0). 



Linearized kinetic equation and moment method 

Frequencies for small amplitude spin-waves in the 
linear regime, such as close to the minimum in 
Fig. 3(a), can be obtained by linearizing the kinetic 
equation For this, we consider small changes 

with respect to the stationary state W^^{z^p)^ namely 
Wmn{z,p,t) = W^^{z,p) ± dWmn{z,p,t). Thus wc in- 
vestigate the spin-waves for short times and small am- 
plitudes. The mean-field, as a function of density like- 
wise expands as Vmn{z^t) = V^^{z) ± 5Vmn{z^t) and we 
substitute both expressions into (|6| without the anti- 
commutator, to obtain the linearized kinetic equation 

{dt-d^)8W =^{[8W,V''] ± [W\3V,]) (8) 

Our next step to explicitly study different kinds of 
spin- wave modes is the so called moment method [43] . 
The thinking behind it is to look at the time evolu- 
tion of moments of the position and momentum op- 
erators. The l-ih. moment of such an operator in 
the phase-space representation is the expectation value 



{z^)mn{^)= J dp J dz6Wmn{z,p^t)zK Expanding the ki- 
netic equation in moments means taking into account 
only different modes of spin-waves up to a certain or- 
der, e.g. dipole (sloshing) modes for / = 1, quadrupole 
(breathing) oscillations for 1 = 2 and further. 

First we need to find a suitable expression for the 
stationary state W^^{z,p), the state of the system 
before a spin- wave is excited by applying a gradient. 
With the preparation scheme in mind, considering the 
preparation pulse to be infinitely short in time, we ap- 
proximate Wo as a product of spin and orbital de- 
grees of freedom W^^{z,p) ^ Mmnfo{z,p) where Mmn 
is a matrix in spin space determined by the prepara- 
tion pulses and /o is the phase-space distribution for a 
two-component Fermi gas onto which the pulse is ap- 
plied. This distribution is given for a non-interacting 
gas in a harmonic trap in local-density approximation 
as fo{z,p)= J d^x J d^p{exp{P{^ ^ lm{ujlx^ ^ uj^y^ ^ 
uj^z'^ — jii)) -\- . In principle, corrections to this dis- 
tribution arise due to interactions but in our case the 
mean-field energy is small compared to kinetic and po- 
tential energy and we neglect them. The mean-field 
V^^{z) = Mmn^oi^) then depends on the density distri- 
bution no{z)= J dpfo{z,p), and we introduce the short 
notation Mmn = ^kii^kinm - Ukmni)Mki- For the fre- 
quencies of the dipole modes we expand 5W and 5V into 
moments of position and momentum up to first order 
only 

6Wnin{z,p,t) = fo{z,p) {Amn{t) ± zBmn{t) ^ pCmn{t)) 

(9) 



SVmn{z,p,t) = no{z) (^Amn{t) zBmn{t) 



(10) 



We take the zeroth and first moments of z^p of 6W in 
([9| and obtain their relationship to A^B^C 

{P)m 



A 
with 



mn — (l) 



mn 1 -Bmn 



{z}rnn 



, Crr 



(11) 



particle number N= J dz J dpfo{z,p)^ 
as well as (^^)o = J dz J dpz^ fo{z^p) and 
{p^)o = / dz J dpp'^ fo{z^p). We substitute expres- 
sions (11) into the linearized equation ([8|, then take 
again the zeroth and first moments of z^p to obtain 
three equations for the matrices A, B, C respectively 



CM 



with 
h = 
h = 



dtB-mujlC=^ ( 

m m 

coefficients 
^ J dz{zno{z))^ 

J dz J dpp^fo{z,p)no{z 



M,A 


± 


A,M 


M,B 


± 


B,M 



(12) 
(13) 
(14) 



lo = J dzno{z)'^ , 
and 

The first equa- 
tion is decoupled from the others and does not lead to 
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spatial dynamics, so we discard it. Equations (13) and 



(14) describe dipole oscillations of all spin components 
of the Wigner function and the frequencies can be 
obtained by a Fourier transform dt — icj and solving 
the eigenvalue equations, similar to the procedure for a 
spin 1 Bose gas in Refs. [34l l35]. 



SU(N) interactions 

For the sake of simplicity we now demonstrate this 
for the special case of SU(N)-symmetry {N = 2s + 1) 
assuming all scattering lengths to be equal. For the 
s = 3/2 subsystem of ^^K considered in the experi- 
ments, our numerical comparison shows only small dif- 
ferences in the spin-wave behavior between this case and 
the true scattering parameters. In the SU(4)-symmetric 
case with all scattering lengths equal, go=g2= • - • =g^ 
the coupling constants are of the particularly simple form 
Uijki = 2 — hi^kj)' For all matrices with a tilde in 

equations ([l2|, ([T3| and ([l4| we obtain the simple ex- 
pression M = ^(Tr(M)l - M). Further, the r. h. s. of 
Eqs. (n!2|) and (p^Sl) are zero and the equations for 5, C, 



written in matrix form reduce to 



dtB - mujlC = 

dtC^-B= ^[M,C] 
m m 



(15) 
(16) 



which decouple trivially. After a Fourier transform dt 
—ioo we obtain the eigenvalue equation 



{u'^ - ujI)C = 2ujuJmf [M, C] 



(17) 



where we have introduced the mean-field frequency 
^mf = ^^2/2^. For all (nematic) initial spin states M con- 
sidered here the solutions of ( p!7[ ) are either the trivial case 
uj= ±ujz that corresponds for instance to oscillations of 
the entire atomic cloud in the trap or 



^ = -^mf ± y CJ^f + . (18) 

which describe the spin-wave propagation. 



S6. TENSOR EXPANSION OF THE WIGNER 
FUNCTION 

The methods described thus far in this supplemental 
material are valid for any value of the spin. We now 
focus on the 5 = 3/2 case in order to demonstrate the self- 
decoupling of the spin- nematic component. In previous 
studies of the s = l/2 case [371 138] the kinetic equation 
and Wigner-function were rewritten in terms of the Pauli 
matrices W = ^{Wot -\- W ■ a) which form a complete 
basis for Hermitian 2 x 2-matrices and allow to identify 
the vector-part of the mean-field as an effective magnetic 



field. For 5 = 3/2 we define a similar basis comprising 
the identity matrix Tq = ^U, the three spin operators 
Tf^ = five quadrupole operators and seven 

octupole operators T^. For the quadrupole we select all 
traceless symmetric matrices quadratic in 5, 









(19a) 








(19b) 








(19c) 


Ti 






(19d) 




"2^3 ^^'^'^ 


SzSy) 


(19e) 



and for the octupole 



41 
20^ 



.^_V5f , 41 



^3 3 

2_V^/^3_41 

^3 - 3 20*^ 

^ 2V^ ^y) 



(20a) 

(20b) 

(20c) 
(20d) 
(20e) 
(20f) 



Tl = {S^SyS, + S.SyS,) . (20g) 

This set of 16 Hermitian matrices forms an ort ho normal 
basis set with respect to the trace 



Tr 



SO we can expand the Wigner function into 

3 21 



(21) 



(22) 



Z=0 m=Q 



with coefficients W{^{z,p) = Ti {Tl^W {z , p)) . This ex- 
pansion is analogous to the moment method shown above, 
but with respect to moments of spin operators. In this 
basis the coupling constants of the mean-field ([t]) have 
the form 



3 21 



cd ' 



(23) 



1=0 m=0 



where the coefficients ai depend on the coupling con- 
stants gs' In a real 5 = 3/2 system, ai =as^ which is a 
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result of the SO (5) symmetry of the system [44]. Analo- 
gous to ([22]) the mean-field can be expanded as 



3 21 



v{z) = Y.^vr{z,p)Tr, 



(24) 



1=0 m=0 



and we find, that in this basis the mean-field is directly 
proportional to the Wigner function with the coefficients 
ai 

Vr{z) = ai J dpWriz^p)- (25) 



We now insert equations (22), (24) and (25) into the 
kinetic equation ([6|, neglecting the anticommutator for 
simplicity and obtain 



3 21 



1=0 m=0 

E / dqWlI^"{z,q)Wir'{z,p) 

l',l"m'm" 



rj-Mm rjim 



(26) 



The time evolution for each tensorial component of 



the Wigner function is obtained from (26) by taking the 
trace with respect to the basis as in (21). On the 
right hand side we define 



^T'm = Tr 



such that S26\j becomes 



rj-im rj-im 



(27) 



{dt-do)wr{z,p) 



ain 



I' ,1" m'm" 



(28) 



Equation (27) simply states the decomposition of the 



commutators of tensor operators in the tensor basis it- 
self. Apart from the trivial case that Tq commutes with 

all other tensors we know that by construction Tl,T( = 

iV^^j^^ijkTi and thus A^-^-^^ = iV^SnCij^- Formulating 
an explicit expression for all values of A[!^//iy"^ is not 
straightforward so we restrict ourselves to provide those 
values of l^V ^l" for which the AJT^^/I^ ^ are non-zero for 
all values of m, m', m": 



1 



(r,F) 



(29) 



(0,0) 

(1.1) , (2,2), (3,3) 

(1.2) , (2,1), (2,3), (3,2) 
(2,2), (3,3), (1,3), (3,1) 



These relations show, which tensor components I" con- 
tribute to the time evolution of the tensor component I. 



To illustrate the these findings in a matrix representation 
we introduce the symbolic notation Wi = Tt{WTi)Ti 
suppressing the m-indices, and insert the relations (29) 
into the equations of motion for Wi. This brings us to 
equation (2) of the manuscript 

dtWo = doWo , 

dtWi ^ doWi + 4 {[WuVi] + [W2, V2] + [W3, Vs]) , 
in 

dtW2 ^ doW2 + 4 i[W2,Vi + V3] + [Wi + W3, V2]) , 
in 

dtWs ^ doWs + 4 {[W3, Vi] + [Wi + W3, V3] + [W2, V2]) , 
in 

(30) 

where the self-decoupling of the nematic spin tensor is 
evident, since no term of the form V2] exists on the 
right hand side of the equation for dtW2. 
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